Semiclassical Theory for Decay and Fragmentation Processes 
in Chaotic Quantum Systems 
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We consider quantum decay and photofragmentation processes in open cliaotic systems in the 
semiclassical limit. We devise a semiclassical approach which allows us to consistently calculate 
quantum corrections to the classical decay to high order in an expansion in the inverse Heisenberg 
time. We present results for systems with and without time reversal symmetry and also for the 
symplectic case, as well as extending recent results to non-localized initial states. We further analyze 
related photodissociation and photoionization phenomena and semiclassically compute cross-section 
correlations, including their Ehrenfest time dependence. 
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I. INTRODUCTION 

Physical phenomena involving decay processes have 
been addressed in many physical contexts. They play a 
central role in the study of excitation relaxation in semi- 
conductor quantum dots and wires [1, 2], photoioniza- 
tion via highly excited atomic [3] or molecular [4] Ryd- 
berg states, photodissociation of molecules [5], atoms in 
optically generated lattices and cavities [6], and optical 
micro-cavities [7], to name a few examples. 

For an open chaotic system it is well known that the 
classical probability of finding a particle inside the sys- 
tem at a certain time, the classical survival probability, 
decays exponentially in time, p'^'(i) = e^*/'^'', where Td 
is the classical life or dwell time. Numerical calculations 
[8], however, revealed that the quantum survival proba- 
bility deviates from the classical one at times compara- 
ble to t* w \/TdtH, where tn = 27rfi,/A is the Heisen- 
berg time, and A is the mean level spacing. Theoretical 
calculations invoking supersymmetry techniques [9, 10] 
confirmed these findings. There it could be shown that 
in the random matrix theory (RMT) limit, the quantum 
decay p(t) takes the form of a universal function, which 
only depends on the general symmetries of the system, 
the classical life time and the Heisenberg time. The first 
successful semiclassical approach to derive the RMT pre- 
dictions for quantum graphs was performed in Ref. [11], 
reproducing the first order RMT quantum corrections for 
networks with and without time-reversal symmetry. 

Recently, we have developed a semiclassical approach 
for calculating the decay of an initially localized wave 
function inside an arbitrary chaotic system [12]. The 
semiclassical framework used there involves correlated 
trajectories which have been shown to be a powerful 
tool and the key to link classical hyperbolic dynamics 
with universal quantum properties [13]. These semiclas- 
sical techniques have been recently extended and widely 
applied in the context of level statistics [14-16], where 
multiple sums over periodic orbits have to be evaluated, 
as well as in the field of ballistic quantum transport in- 



volving Landauer-Biittiker formulas [17-23], where tra- 
jectories start and end at the openings where the chaotic 
conductor is attached to leads. In Ref. [12] a unitarity 
problem was encountered when using these semiclassical 
techniques to evaluate the contribution of pairs of inter- 
fering trajectories starting and ending inside the system. 
Therefore a new kind of diagram was considered, which is 
crucial for ensuring unitarity in problems involving open 
trajectories connecting two arbitrary points in the bulk. 
A similar type of trajectory appears in the semiclassi- 
cal description of transport if the coupling between the 
chaotic conductor and the leads is not perfect, as shown 
in Ref. [24]. 

In this article we generalize the approach presented 
in Ref. [12] for localized initial wave functions to non- 
localized wave functions. We outline how to systemati- 
cally obtain higher order (in t/t*) quantum corrections 
to the classical decay and present terms up to the 7th 
order and 8th order, for systems with and without time 
reversal symmetry, respectively. We further calculate the 
survival probability for systems with spin-orbit interac- 
tion, corresponding to the symplectic RMT ensemble. 

Closely related to quantum decay are problems of 
atomic photoionization or molecular photodissociation 
where the fragmentation mechanism involves photoex- 
citation to an intermediate excited resonant state (with 
corresponding complex classical dynamics) which then 
subsequently decays by sending out a particle, i.e. an 
electron, atom or an ion. In the semiclassical limit, spec- 
tral correlation functions for the related photoionization 
and photodissociation cross-sections can be expressed 
through the spectral form factor and the survival proba- 
bility. Earlier semiclassical treatments [25, 26] of photo 
cross-sections were always limited through the diagonal 
approximation used which was relaxed in this context 
only very recently [12]. Here we will present a detailed 
semiclassical treatment of the brief account given in [12] 
and extent the results by including Ehrenfest time effects 
for cross-section correlations and by computing higher or- 
der contributions. 
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This article is organized as follows: in Sees. II and III 
we present the semiclassical approach to the quantum 
survival probability, generalized to non-localized wave 
functions, by including a time average. In Sec. IV this 
approach is further extended to derive higher order cor- 
rections for systems with and without time reversal sym- 
metry as well as for the case of spin-orbit interaction 
which follows the universal RMT prediction for the sym- 
plectic case. In Sec. V we analyze fluctuations of the sur- 
vival probability through its variance. In Sects. VI and 
VII we give a detailed semiclassical analysis of the statis- 
tics of photofragmentation, including higher order correc- 
tions and the Ehrenfest time dependence of the leading 
quantum contributions. We conclude with an outlook in 
Sec. VIII. 



II. SEMICLASSICAL APPROACH TO THE 
SURVIVAL PROBABILITY 

The quantum mechanical survival probability as a 
measure of the decay is defined as 

p{t)= [ dri,{r,t)r{r,t), (1) 

J A 

where ^(r, t) is a wave function and A the volume of the 
system we are considering. For a closed system p{t) = 
1, while for an open system this no longer holds and 
p{t) decays in time. Expressing ip{v,t) in terms of the 
propagator K{r, r', t), 



i>{v,t)= / dv'K{v,v',t)M^')., 



(2) 



we have 



p{t)= f dvdv'dv"K{v, r', t)K*{v, r", t)Vo(r')V'o(r") , 

J A 

(3) 

where "00 (r) is the initial wave fimction at f = 0. 

In order to calculate the semiclassical expression for 
p {t), we replace the exact quantum propagator K{r, r', t) 
with the semiclassical Van Vleck propagator [27], 



K'''' (r,r',t) 



1 



^ ' 7(r'^r,f) 



(4) 



Here / is the dimension of the system (in the following 
we will consider / = 2), S^(v, r',t) = dt' Lj[tj,rj, t'] is 
the classical Hamilton's principal function (with the 
Lagrangian) along the path 7 connecting r' and r in a 



time t, and Dx, 



det 



dvdr' 



1/2 



e '2^*7 is the 



The semiclassical survival probability is then given by 

= 7^ / dvdv'dvMr'Woir") X (5) 

X ^ D;^Li|,e^('S^-Sy). 

7(r'— 'r,t) 

In the following, we introduce a local time average 
in the survival probability which enables us to neglect 
highly oscillating terms in the above double sum. We 
deflne 



Pit) ^ (p^^(i))At ^ — / P'''{t')dt' 



(6) 



with At ^ t. We will later see that for a localized ini- 
tial wave packet p{t) w p{t) in the semiclassical limit, 
recalling the result of Ref. [12]. 

The phase difference in the double sum in Eq. (5) 
rapidly oscillates unless the two related trajectories are 
correlated. Therefore most of the contributions will dis- 
appear due to the time average. The contributions that 
prevail the average are from pairs of correlated trajec- 
tories with action differences of the order of h, which 
implies that the trajectories 7 and 7' should be 'similar'. 
This puts a restriction on the initial points of the two 
trajectories, i.e. they should be almost the same. We can 
then expand trajectories 7 (or 7') going from r' (or r") 
to r in a time t around trajectories 7 (or 7') going from 
ro = (r'-|-r")/2 to r in a time t. This expansion amounts 
to approximating the classical prefactors D^(r, r',f) « 
D^{Y,YQ,t) and £'^'(r, r",t) « D-^/ (r, ro, t), while ex- 
panding the phases in the exponents up to the first order, 
because the latter are more sensitive to small changes in 
their argument. The expansion of the actions yields 



(r, r', t) « (r, tq, - • Pt.o , 



(7) 



S^, (r, r", t) « S^, (r, tq, + -q • p^,,o, (8) 



where q = r' — r" and p-y_o (or P7',o) is the initial mo- 
mentum of the trajectory 7 (or 7'). The semiclassical 
survival probability, Eq. (5) then reads 

= ( (d^ / '^'•^'•o'^q^o + f ) V'o (ro - f ) X 

X D^£);,e*(^^-^v)e-*P°V'>) , (9) 

7,7'(ro— >r,t) 

where , = (p7,o + P7',o)/2. This can be written as 



Pit) 



1 



drdra x 



Van Vleck determinant including the Morse index /ij. 



(27rfi)2 

X J2 D,D;,e^^^^-^y^pw {ro,P°^,) 

7,7'(ro^r,f) 



(10) 



At 
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where 



Pw(r,p) 



r + Y ) V-o 



.r p 



(11) 

is the Wigner transformation of ipoir). For an initial 
coherent state, the integrals over Tq and r' can easily be 
performed, and the result is consistent with that of Ref. 
[28]. 

Eq. (9) still involves rapidly oscillating phases, and 
again most of the contributions will cancel out, unless the 
trajectories in a pair are systematically correlated. The 
main contribution corresponds to the diagonal approx- 
imation, i.e. 7 = 7', which gives the classical survival 
probability. Together with the sum rule [29] for open 
systems, this yields 



P 



where 



p indicates a phase space average, 



J drdpF{r,p)pw{r,p) 



(12) 



(13) 



and 1/tij is the classical escape rate at the energy E — 
H{r, p), where H{r, p) is the Hamiltonian of the system. 
For a two-dimensional system, ~ fl{E)/{2wp), with 
n{E) = J dv'dp'5{E - H{v', p')), w the size of the open- 
ing, and p — |p|. For a chaotic billiard this reduces to 
Td = rmrA/(wp). For an initial state with a well defined 
energy Eq we can write p'^''''^{t) = e-*/^d{Eo)^ ^^ic fol- 
lowing we will assume this to be the case and drop the 
brackets of the phase space average. 

Equation (12) has two restrictions: First, we have sup- 
posed that at time t the trajectories can already be con- 
sidered ergodic (they have homogeneously explored the 
phase space). This is a good assumption as long as 
tX ^ 1, with A being the Lyapunov exponent. Second, 
we have assumed that the ergodicity of the corresponding 
closed system is not affected by the opening, meaning, 
classically the opening should be small r^A 3> 1, while 
quantum mechanically it is very large ^ tn- 



III. SURVIVAL PROBABILITY: LEADING 
ORDER WEAK LOCALIZATION-TYPE 
CONTRIBUTIONS 

It was shown in Ref. [12] that the leading quantum 
corrections to the semiclassical survival probability (5) 
for systems with time reversal symmetry come from or- 
bits with a self encounter (Fig. la), 'two-leg-loops' (211 or 
2-encounter) introduced in Ref. [13], together with fone- 
leg-loops' (111, sketched in Fig. lb, c), which together 
preserve unitarity. 
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FIG. 1: (color online) Scheme of two-leg-loop (211, a) and one- 
leg-loop (111, b,c) orbit pairs. The trajectories 7 (full line) and 
7' (dashed) connect the points ro with r in a time t, and they 
differ by a 2-encounter in (a). When the beginning or the end 
of the trajectory is inside the encounter we have the situation 
plotted in (b). (c) is a variation of (b) where there is no 
self-crossing of either of the two trajectories. 



Two- leg-loops 

In this section we will give a detailed derivation of these 
contributions to the survival probability following the 
phase space approach [15]. The double sum over trajec- 
tories is replaced by the sum rule together with integrals 
over the stable and unstable manifolds along reference 
trajectories 7 weighted by the density of 2-encounters in 
a orbit of length w^''(m, s,t), giving rise to a difference 
in action AS'(m, s) = ms, whose absolute value is smaller 
than a classical value (? . This density is given by 



(t - 2ic 



2^^ 



(14) 



where the encounter time is icnc — A^^ ln(c^/|Ms|). 

The classical survival probability is modified by a fac- 
tor e*""'^/'^'', since the fact that the first stretch remains 
inside the cavity implies that the second will also be in- 
side. Thus 



du\ dsw2ii(y^s^i)e*™c/r,gt«.^ (^5^ 



The integration can be performed by making the change 



of variables 
result is 



— usjc^, cr = c/u as in Ref. [20]. The 



p--(t) 



e t 
2— 



(16) 



The quadratic term corresponds to the first order quan- 
tum correction according to Ref. [9] , while the linear term 
breaks unitarity, since it does not vanish as — > 00 
(when the system is closed). As shown in Ref. [12] an- 
other type of diagram has to be considered in order to 
solve this problem. 
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One-leg- loops 

The relevant diagrams correspond to trajectories with 
an encounter at the beginning or at the end of the tra- 
jectory, as shown in Fig. lb, c. Clearly, the latter only 
exists for initial and final points inside the cavity, since 
at the openings the exit of one stretch of the encounter 
implies the exit of the other one (with perfect coupling) . 

To evaluate these two contributions we define a 
Poincare surface of section at some time t' from the end 
or beginning of the trajectory [20]. The encounter time 
will be given by 



t,nc{t',u)^t' + j\n{c/\u\).. 



(17) 



with the restriction t' < j\n{c/\s\), while the density of 
such encounters is given by 



w^^\u,s,t) = 2 / dt' 

Jo Jo 

■ In T% 



dh 



dt 







1 1 2icnc (^ 1 ^) 



(18) 



The factor two is due to the possibility of having the en- 
counter at the beginning of the trajectory or at the end. 
The difference in action will be AS ~ us at any point of 
the Poincare surface of section. It is important to men- 
tion that this weight function automatically includes the 
situation where both end points are very close i.e. co- 
herent back-scattering. We can now proceed to calculate 
this contribution to the survival probability in the same 
way as before, replacing w^''(m, s, t) by w^^^{u, s, t) in Eq. 
(15). In order to evaluate the integrals, we make the 
change of variables [20] 



= c/cr, s = cxa, (19) 



with an integration domain — l<a;<l,l<o'<e 
and < t" < y In . Here is important to notice 

that the limits of t" also include the situation where the 
point at which the orbits start is after a possible self- 
crossing. This means that it is not necessary to have a 
true self-crossing in configuration space in order to give 
a contribution of this kind. 

We define p^^^{t) = /e^*/'^'' where. 



7 = 2/ du I ds 

c J —c Jo 



i — Of 

^^cnc 



(20) 

the integral over a can be easily done after the change of 
variables mentioned above, and / can be written as 



4rA 
T^tn Jo 



1 /-ilnll/x) 

dxcos{rx) / dt"{t-2t")e^ 







4rATd 
drT^ J irtn 



dxcos{rx)x 



(21) 





FIG. 2: (color online) Scheme of orbit pairs that do not re- 
quire time-reversal symmetry that give higher order correc- 
tions: (a) a single 3-encounter, (b) a double 2-encounter. The 
trajectories 7 (full line) and 7' (dashed) connect the points 
ro with r in a time t, and they differ by the way they are 
connected in the encounter regions. 



where r — (? jTi. 

The integration over x can be performed by parts, ne- 
glecting highly oscillating terms that will disappear after 
averaging [20], yielding 





2 ' 
















f dy 

JO 



Ar\Td ( sin(r) 



TTtH 

sin(j/) _ 

y ~ 



it 
T^tH Jo 



dy 



1 

sin(y) 



dx 



sin(ra;) 



y 



(22) 



I'M 



Then the 111 contribution to the decay reads 

t 



p-m it) = 2-e 

tH 



(23) 



This term exactly cancels the linear term in Eq. (16) 
coming from the 211 contribution, recovering unitarity. 
The leading semiclassical correction (quadratic in time) 
to the classical survival probability is therefore [12] 



P' 



■211+111 



2Tdt 



H 



(24) 



which is consistent with the RMT prediction [9]. It can 
be interpreted as a interference-based weak localization- 
type enhancement of the survival probability. 

In the next section we will extend this approach to in- 
clude higher order corrections, coming from semiclassical 
diagrams with multiple encounters or with one encounter 
involving multiple stretches. 



IV. SURVIVAL PROBABILITY: HIGHER ORDER 
CONTRIBUTIONS FOR THE GUE, GOE AND 
GSE CASES 



For the unitary case, the next order contributions to 
p{t) are given by the diagrams shown in Fig. 2, as in- 
dicated in Ref. [15]. In a similar way, we can compute 
the next order corrections for systems with time rever- 
sal symmetry. Time reversal symmetry, however, allows 
more structures, the corresponding diagrams include the 
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ones sketched in Fig. 2 (multiplied by a factor of 4 for 
Fig. 2a and a factor of 3 for Fig. 2b, respectively [15]) 
together with a structure including two copies of the en- 
counter in Fig. la. 

In general, an encounter region contains an arbitrary 
number of Z > 2 stretches of the trajectory, which are 
mutually lincarizablc, and one speaks of an Z-cncounter. 
In order to calculate higher order corrections, we consider 
trajectory pairs with encounters described by the vector 
V, whose elements vi list the number of Z-cncountcrs in 
the trajectory pair. The total number of encounters is 
then V = ^vi while the number of links of the related 
closed orbit is L = ^ Ivi as in Ref. [19]. 

Consider a periodic orbit formed by joining the ends 
of the open orbit. We can generate the open trajectories 
by cutting this closed orbit along each of its links and 
moving the ends of the cut to the required positions. 
Note that for systems with time reversal symmetry, we 
must chose cither the partner orbit or its time reversal 
so that the link, which is cut, is traversed in the same 
direction by both orbits. The contribution can then be 
separated into three parts: 

A where the start and end points are outside of the en- 
counters (211), 

B where either the start or end point is inside an en- 
counter (111) and 

C where both the start and end point are inside encoun- 
ters (Oil). 



Case A 



This contribution can be written as 

Pv,A{t)=N{v) j As dti w.v,A{u,s,t)t 



-jit 



while the last link time is fixed by the total trajectory 
time 



L+l V 
i=\ a=l 



(27) 



When we perform the integrals the weight function be- 
comes 



111 A (ll t t) ~ ft Igtcnc) 



(28) 



To calculate the semiclassical contribution we will rewrite 
Eq. (25) as 

Pv,A{t) = ^{v) j ds du z„,a{u, s, ^)e-''*e^"^ (29) 

where Zv,a{u, s, t) is an augmented weight including the 
term from the survival probability correction of the en- 
counters 

z^^Aiu, s, t) = w^^a{u, s, t)e^»('"-i)''*?„c (30) 

where we have expanded in the second line the exponent 
to first order in the encounter times. We can now use 
the fact that the semiclassical contribution comes from 
terms where the encounter times in the numerator can- 
cel those in the denominator exactly [15]. Keeping only 
those terms, we then obtain a factor of {2-Kh)^~^ from 
the integrals over s and u and obtain the result for tra- 
jectories described by the vector v of interest. 

Consider for example a trajectory with a 3-encounter 
with two long legs, sketched in Fig. 2a. The encounter 
has a duration given by 



X e 



Er=i('c-i);it?„cgiMs 



(25) 



In- 



(31) 



where N(y) is the number of trajectory structures cor- 
responding to each vector v [15], /z = l/r^, and a labels 
the V encounters, each being an la- encounter. We have 
included the correction to the survival probability of the 
trajectories due to the proximity of encounter stretches 
during the encounters. In terms of an integral the weight 
is given by 



Wv,a{u, S, t) = 



St 



dt. 



■Jo 



t-tcnc 



-f2 



dti 



(26) 



where tcnc is the total time that the trajectory spends 
in the encounters t^nc = ^^=1 ^a^nc- Each of the links 
must have positive duration and this restriction is in- 
cluded in the limits of integration. The weight is simply 
an L-fold integral over different link times ti,i = 1 . . . L, 



where j = 1,2 and the differences between the 

unstable and stable coordinates of the trajectory on the 
PSS placed in the encounter region, respectively. 

The density of this type of encounter, with an action 
difference AS = u ■ s, is 



(t 3ienc)'^ 
enc 



(32) 



where we use the notation (I)"' to indicate that the tra- 
jectory has vi Z-encounters. We can calculate the contri- 
bution of such orbits by replacing the sum over the part- 
ner trajectory 7' with an integral over the stable and un- 
stable coordinates (tt, s) with the density «^(3)i,a('W, s, t), 
modifying the classical survival probability entering the 
sum rule by a factor e^'^*""" . In the case of time reversal 
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symmetry there are four possible structures in this case 
[19], and the final result is 



P(3)^A(^)=4e-*/-''(-:^ 



(33) 



For a double 2-encounter shown in Fig. 2b, we define two 
encounter times: tl^^ w i In , , and t^„^ w i In , , . 
The density of such a double-encounter is given by 



{t - 2te 



(34) 



with tcnc = ^onc+*cnc- ^^lis casc the number of possible 
structures for systems with time reversal symmetry is 5. 
The contribution of such orbits to the survival probability 
is 



P(2)^A(^) = 56-*/^'' (2 



2*3 



■ (35) 



The total contribution of structures with L — V = 2 of 
2U's is then 



Pi 



,A(t) = e-*/^-* (4 



2e Si^ 
+ 



(36) 



Case B 

Now we have to consider the corresponding 'one-leg- 
loops' for the previous diagrams. This contribution can 
be written as 

PvA^) = ^(v) j ds du ^„,b(w, s, t)e-''*es"^ (37) 

Here one encounter overlaps with the start or end of the 
trajectory, we have therefore one link fewer (L in total) 
and an extra integral over the position of the encounter 
relative to the starting point. Starting with a closed peri- 
odic orbit, (and dividing by the overcounting factor of L) 
we can cut each of the L links in turn and move the en- 
counter on either side of the cut to either the start or the 
end. In total we obtain copies of the same 111 involv- 
ing the encounter a', and additional factor of 2 appears 
due to the possibilities of having the encounter at the be- 
ginning or at the end of the trajectory. The augmented 
weight can then be expressed as a sum over the different 
possibilities, each of which involves an integral over the 
distance from the PSS to the initial or final point, *„/, 



L-l 



xe^«=i('o-i)''*?nc. (38) 

Because of the integrals over the position of the encounter 

at the start or end of the trajectory, the semiclassical con- 
tribution is calculated differently, using integrals of the 



type we encountered in Eq. (21). However, it is easy 
to see in Eq. (21) that after a suitable change of vari- 
ables, the integral over r can be effectively replaced by 
a tcnc- This change of variables can be done for each 
(wq,/, Sc', ia')> giving again a factor of f^^^ for each inte- 
gral over ta' , so that the augmented weight can be written 
as 



Zv,-B{u,S,t) ' 



, ^ (X^g ^g^cnc) X^ct ^«^enc) 

n(i+(^a-ix„c), 



(39) 



and treated as before. 

For a single 3-encounter, we define again a Poincarc 
surface of section at a time t' from the beginning or end of 
the orbit. The encounter time is given by tenc(i', Wmax) = 
t' + \ ln(c/maxi |ui|), and the augmented weight of such 
encounter 



f I s Imax 

/ dt 
Jo 



(40) 



Making the change of variables indicated in Eq. (19) and 

multiplying by the number of possible structures, the re- 
sulting contribution for systems with time reversal sym- 
metry is 



P(3)i,B(i)=4e-*/^'^(|;). 



(41) 



The integration of Eq. (40) yields the same result as 
if we had used instead the augmented weight function 
(t-3«enc)^ (^ ^ 2^te„c), as in Eq. (39). 

For the double 2-encounter one of the Poincare surfaces 
of sections will be from the beginning (or end) of the 
trajectory at a time t', then tl^^i^' , Ui) =t' + jln \c/ui\, 
and the corresponding density of such pairs is 



1 



■ In - 



dt 



2*cnc)-'^^t_^ (42) 



''enc''enc 



We perform the same change of variables as before for 
(ui, Si, r), which yields 



P(2)2,B(t) = 5e-*/"'' 



'MjjTd t 



H 



(43) 



The total contribution of Ill's iov L—V = 2 for systems 
with time reversal symmetry is given by 



P2,B(i) 



(44) 



Case C 

This contribution can be written as 

Pv,c{t) = j dsduz^,c{u,s,t)e-^'*e^^\ (45) 
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Now that we have one encounter overlapping with the 
start of the trajectory, and a second (different) encounter 
at the end of the trajectory, we have several additional 
complications. Firstly, there is again one link fewer {L— 1 
in total) and now we have two extra integrals over the 
position of the start and end encounters relative to the 
start and end point. Also the number of such structures 
is different. Starting with a closed periodic orbit, we can 
cut each of the L links in turn and move the encounters 
on either side of the cut to both the start and the end, as 
long as the link joins two different encounters. We there- 
fore need to count the number of ways that this is possible 
for the different sizes of encounters that are linked. We 
record these numbers in a matrix J\f{v), where the ele- 
ments Afa,0{v) record the number of links (divided by L) 
linking encounter a with encounter (3, in this case it is 
convenient to include J\fa,/3{v) in the augmented weight 
function. The augmented weight, including these possi- 
bilities, can then be expressed as the following sum over 
the Oil encounters 



Zv,c{u,S,t)=^J\fa',f3'{v) / dta'dtfj' 
a',l3' ■' 



.Er=i('c-l)Mte° 



(46) 



(L-2)!0^-vn, 

Again we can expand the exponent to first order in the 
encounter times and write the augmented weight function 
as 



1^2 



(L-2)!Q^-^nai?nc' 
and treating it as before. 



(47) 



For a single 3-encounter there cannot be such a con- 
tribution. For a double 2-encounter we may define two 
Poincare surfaces of section at t[ and from the be- 
ginning and the end of the trajectory, ienc(*'i' — 
t[ + iln|c/ui| and tonc(*2)'"2) = t2 + xlii|c/w2|, and 
the corresponding density 



■ In - 



(t 2tenc) 



(48) 



This gives the following contribution to the survival prob- 
ability 



P(2)2,c(0 = e 



(49) 



Again the result is the 
to Eq. (47). 



same 



2^2 

''H 



as if we had used 
nc), which corresponds 



UnitEiry case 



We can easily calculate the contribution for each vector 
V for each of the three cases, as long as we know the num- 
bers of possible trajectory structures. For cases A and B, 
these numbers can be found in Ref. [15] and are repeated 
in the first four columns of Table I. For case C we will go 
up to the sixth order correction, L — V = 6, and for this 
we have at most three different types of Z-encounters. It is 
useful to rewrite the sum over a and /3 as a sum over the 
components of the vector v. Afa,fj{v) records the number 
of ways of cutting links that connect encounter a and 0, 
in the periodic orbit structures described by v. However 
we can see that the important quantities are the sizes of 
the encounter a and /3. Instead we record in Nk,i{v) the 
number of links that join an encounter of size k to an 
encounter of size I. If we number the encounters from 1 
to V in order of their size, then we only need to know 
the numbers Ni^^i^(v), Ni^^iy{v) and Niv_^^iy{v), as the 
maximal number of different sized encounters is three. 
Moreover J\fk,i is symmetric, therefore we include in Ta- 
ble I both A4,i and Ni^k together. Using a program to 
count and classify the possible permutation matrices we 
obtain the remaining columns in Table I for systems with- 
out time reversal symmetry. Note that certain encounter 
combinations might correspond to several elements of the; 
numbers Muhi^), ■N'h,iv{v) and Miy_^,iy{v), in which 
case we record their number in the leftmost column. 



V 


L 


V 


N{v) 








{2Y 
(3)1 


4 
3 


2 
1 


1 
1 


1 






(2) - 
(2)2(3)1 
(2)1(4)1 

(3) 2 
(5)^ 


8 
7 
6 
6 
5 


4 
3 
2 
2 
1 


21 

49 
24 
12 

8 


21 
12 
16 

8 


32 




(2)° 
(2)*(3)i 


12 


6 


1485 


1485 






11 


5 


5445 


2664 


2592 




(2)^(4)1 
(2)2(3)2 
(2)2(5)1 
(2)1(3)1(4)1 


10 


4 


3240 


984 


1920 




10 


4 


4440 


464 


2624 


960 


9 


3 


1728 


228 


1080 




9 


3 


2952 


552 


760 


1080 


(3) ^^ 

(2) 1(6)1 

(3) 1(5)1 

(4) 2 
(7)1 


9 


3 


464 


380 






8 
8 
8 


2 
2 
2 


720 
608 
276 


360 
360 
180 






7 


1 


180 









TABLE I: The number of trajectory pairs and the number 
linking certain encounters for systems without time reversal 
symmetry. 



Table I allows us to obtain the following results for the 
quantum corrections to the classical decay for the unitary 
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case 



P2{i) 
Pi{t) 



7- / +4 



t 



H 



ISOrJ 



1920t4 



t% V224rJ 



+ 



31t 



10 



22680rf 



30240r^ 



10080Tj 322560t5 



(50) 
(51) 

(52) 



There the symplectic case is obtained by inchiding in the 
Hamihonian a classically weak spin-orbit interaction. 

In the following we study the effect of spin-orbit inter- 
action on the survival probability. The spin-orbit inter- 
action is accounted for by replacing the Hamiltonian for 
the orbital dynamics, Hq considered up to now, by 



fl' = Fo + s-C(x,p), 



(55) 



These results enable us to calculate the decay up to 
8th order in t, giving as the final result 



with C (x, p) characterizing the coupling of the transla- 
tional degrees of freedom to the spin operator s. 



-GUE 



1 + 



+ 



1920Ti4 



Orthogonal case 



180t|4 



For weak spin-orbit interaction, the semiclassic;al prop- 
agator is similar to Eq. (4), where the classical trajecto- 
ries are the same as for the case without interaction [33] . 
. (53) The only modification appears in the prefactor that 
contains now the additional factor (x', p', t), which is 
the spin-,s representation of the spin propagator matrix 
(x'j p', t), defined as the solution of [33], 



Similarly, we can find all possible permutation matrices 
and obtain Table II (sec Appendix A) for systems with 
time reversal symmetry. This gives us the result up to 
7th order in t 



-GOE 



it) 



I 



'i.TdtH 



(54) 



H 



+ 



I 



11 



+ 



/ 41 



2 

7 



V720riti, 
29 
168^ 



15Td4 



+ 



14 



+ 



16 



2lTdt% 



t' +... 



The predictions for the decay using supersymmetry tech- 
niques can be found in Rcf. [10], where the integrals ap- 
pearing there can be expanded in powers oit/tn, follow- 
ing the steps indicated in Ref. [30]. The results of these 
expansions agree with Eqs. (53) and (54). 



d 



^ (x', p', t) + -a-C (X(t), P(t)) 6^ (x', p', i) = , 



with the initial condition 6^(x',p',0) = 1. 



(56) 
This prop- 



agator can be used now in the derivation of a modified 
formula for the survival probability in the case of spin 
orbit interaction. After replacing the initial state [V'o) 
introduced in Eq. (2) by [*o) = \i^o® sq) containing ad- 
ditionally the initial spin state |so), we obtain the matrix 
element (sq |i?-yi?^| sq) as an additional factor inside the 
double sum in Eq. (9). We are interested in the aver- 
age behaviour of this quantity. Therefore, we analyse 
■^25q3YyTr(i?^i3t ) with Tr denoting the trace in the spin 

space. This quantity was already considered in Ref. [32] , 
where by assuming the mixing property of the combined 
spin and orbital dynamics, i.e. full spin relaxation, it was 
shown that we can effectively write 



(-1) 



2s 



L-V 



2s + l 



(57) 



Spin-orbit interaction and the symplectic case 



Along with the cases with and without time reversal 

symmetry, there has recently been interest in a semiclas- 
sical treatment corresponding to the symplectic RMT en- 
semble in different contexts, such as in spectral statistics 
[15] and in the quantum transmission through mcsoscopic 
conductors in the Landauer-Biittiker approach [31, 32]. 



with L and V defined as before. It is important to notice 

that the contribution from spin orbit interaction depends, 
apart from on the spin quantum number s, only on the 
difference L — V. The term (57) can now be inserted 
as prefactor into the expressions in appendix A for the 
GOE case after choosing in each term the correct value 
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oiL-V. For s = 1/2 this yields 



function of the initial state, thus 



-GSE 



it) 



4:TdtH 



(58) 



+ 



H 



H 

A 



( 11 



H 



41 



29 

2688r|t|^ 



192x^4 



240t24 



+ 



84Td4 



This result is again consistent with RMT-type results 

for the symplcctic ensemble [9]. The second, negative 
term in (58) reflects weak-antilocalization effects in the 
quantum decay. 



V. VARIANCE OF THE DECAY 

In Sec. 11 we introduced a local time average, in order 
to select from the trajectories contributing to the double 
sum in Eq. (5) those that start from the same point. 
In order to compare deviations of pif) from the time- 
averaged /5(f), we consider on the level of the diagonal 
approximation, the variance of p(t), averaged again over 
a time window: 



(59) 



Substituting Eq. (4) in Eq. (59), we can write the vari- 
ance as 



(ri)V'o(r2)V'o(r4)V'o(i"5) 



E 



7l(ri^r3,t) 73(r4->rg,t) 
T'2("'2-»''3>*) ')'4("'5-*'"6>*) 



^e^ 



At 



(60) 



where A = D^y^Di^D^^Di^ and AS = S^y, - Sj^ + - 
S^^. Here the configurations ri w r2 and r4 « have 
already been taken into account in p{t)^ and therefore 
have to be ignored in Eq. (60). Due to the average most 
of the contributions to Eq. (60) will cancel out. so for 
surviving the average the configuration of the points 
must be such that the phase difference AS tends to zero. 
Apart from the configurations that contribute to p(t) , the 
leading contribution comes from ri « and r2 « r^, 
which requires w re. We expand the trajectories 71 
and 74 around trajectories 71 and 74 going from qi = 
(ri -|-r5)/2 to qs = (r3-|-r6)/2 and trajectories 72 and 73 
around trajectories 72 and 73 going from q2 = (r2-|-r4)/2 
to qs- We can perform the integrals over ri — and 
r2 — r4 and write the variance in terms of the Wigner 



Sp\t) = 



{2'KhY 



(61) 



!=1 



')'1.74("ll->13.*) 
')'2.73 (12-* 13.*) 



with 



= _ (Pti.O + P74,o) -0 

P7I74 Q ) P7273 



(P72,0 + P73,o) 



and 



Here 



5S. 



, (62) 



(63) 



= (P71,/ - P72,/ - P73,/ + P74,/) • (rs - re)/2 (64) 

where P7i,/ stands for the final momentum of trajectory 
7i- 

We consider here only the contribution from the diag- 
onal terms 71 = 74 and 72 =73, which leads to 



Sp\t) 



1 



diag ■ 



i=l 



E 1^7.1 

7i(qi-^q3'*J 



xen^Sip^ (qi,p^^,o) PW (q2, P72,o) )^^- (65) 

Here A5d = (P71,/ - P72,/) • (rs - re). Upon applying 
the sum rule [34] this can be written as 



5p\t), 



1 



{2-KhYA 



dk 



Al 

(66) 

For a Gaussian initial state, the integrals can easily be 
performed, and for Xt ^ 1 we obtain 



2'Kah 
Apo 



0, 



(67) 



where a denotes the spatial width of the initial state and 
Po the magnitude of its mean initial momentum. Here a 
few remarks about are due: (i) The result in (67) should 
be considered as an estimate of the leading-order h con- 
tribution to the variance as it is based on the diagonal 
approximation. The fact that it is not strictly zero in the 
limit ^ 00 (closed system) makes \is believe that there 
are further contributions, cancelling this term for the 
closed system, (ii) Equation (67) describes 'mesoscopic' 
fluctuations of the survival probability which turn out to 
be non- universal as (5/0^(f)diag scales with the width a of 
the initial state [35]. (iii) Expression (67) may explain 
decay fluctuations which have been found from numeri- 
cal calculations of the quantum decay based wave packet 
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propagation [12]. (iv) Furthermore, for a localized wave 
packet 6p^{t) ^ as ^ 0, and we have p{t) ~ p{t), 
recalling the result in Ref. [12]. 

(v) The variance (67) can alternatively be written as 
(5/0^(<)diag ~ e~^*/'^''/M, where M is the number of eigen- 
states of the closed system necessary to expand the initial 
wave function. 



VI. STATISTICS OF PHOTOFRAGMENTATION 
CROSS-SECTIONS 

Typical examples of quantum decay processes are 
molecular photodissociation [4, 5] or atomic photoioniza- 
tion [3, 36], where the molecule (or atom) absorbs one or 
several photons such that the system is (highly) excited 
to an intermediate configuration coupled to the contin- 
uum, which subsequently allows for decay, i.e. dissocia- 
tion or ionization of the system. 

If this decay is sufficiently slow, a large portion of the 
complex, presumably chaotic phase space of the excited 
system can be explored and the statistics of such pro- 
cesses are assumed to show universal behaviour, as de- 
scribed by the RMT approach developed in Refs. [37, 38] . 
In these "indirect processes" the effective Hamiltonian of 
the excited molecule or atom can be written as H — iT/2, 
where H is the Hamiltonian that represents the part of 
the Hamiltonian containing the "binding" potential and 
r is a matrix describing the coupling of the system to N 
external open channels, which are the possible states of 
the dissociated molecule (or remaining ion). 

Semiclassical approaches to the auto-correlation func- 
tion of photodissociation cross sections were still limited 
by the diagonal approximation used in Refs. [25, 26], 
which however adequately describes the leading order in 
1/N Lorentzian profile of the correlation function. In 
Ref. [12] we briefly presented the leading off-diagonal 
quantum corrections for systems with time reversal sym- 
metry. The purpose of this section is to develop a 
semiclassical approach for quantum corrections to the 
photofragmentation cross-section for systems with and 
without time reversal symmetry, including higher order 
corrections and finite Ehrenfest time effects. We follow 
the diagrammatic approach in Ref. [25] and introduce 211 
and 111 contributions in order to calculate the quantum 
corrections. We will see that the form factor of the cross- 
section auto-correlation function can be semiclassically 
written as the the sum of the survival probability based 
on open trajectories in the excited system (and weighted 
by a factor which accounts for the symmetry) and the 
spectral form factor related to periodic orbits that remain 
trapped inside the system. 

We consider the disintegration of a molecule from its 
ground state \g) via photoexcitation through an interme- 
diate excited electronic surface. 

The photodissociation cross-section of the molecule, in 



the dipole approximation, is given by [5] 

a{E) = ImTr{AG- (E)} (68) 
= ImjdrJ dr'A{r,r')G-{r',r,E), 

where G~{E) is the retarded Green function of the 
molecule, A is a projection operator, given by 

i = 77|</>)(0|, \4')=D\g), (69) 

where = d • e is the projection of the electric dipole 
operator of the molecule, d, on the polarization axis e of 
the absorbed light, and ry = E/chtQ. 

The two-point correlator of the cross-section is defined 

as 

, _ {a[E + h^/2)a{E - h^/T)) - {a{E))^ 
C{u;)= ^^^^p , (70) 

where (...) denotes a local average in energy around E 
and {(j{E)) is the mean cross-section. In the semiclassical 
limit {a{E)) w a{E), where 



a{E) 



(27rfi)2 



drdpAw{r,p)6{E-H{r,p)), (71) 



with 



,p)= j dr' 



Aw{r,p) = / (ir'(r + r72|i]r - r72)e-^^ 



(72) 



the Weyl representation of the operator A. 

In the following, we consider the Fourier transform of 
C{lu), the cross-section form factor. 



Z{t) [ dwe'^'Ciu). 

27r i_oo 



(73) 



As C{ui) — C{—Lo) then Z(t) is real and even. We con- 
sider Z{i) for t > and calculate C(a;) from C(a;) = 
^/o°°Z(i)cos(a;t)rft. 

In order to calculate the semiclassical expression for 
this quantity, we replace the exact Green function by its 
semiclassical counterpart [27] , given by 

G^^ {v\v,E) = ^Te^^^^"'"" (74) 



j{r—>r',E) 



for a two-dimensional system, where D^^ = 



det I — J, 1 exp (— and is the 

Morse index plus additional phases (see Ref. [27]) and 
5-y(r, r', i?) = p.y • dq^ is the action integral along the 
trajectory 7 connecting the points r' and r with fixed 
energy E. 

The semiclassical cross-section form factor (73) is then 
given by 

^"^^^ = 8;?fc^'(/ n'^rM(ri,r2)A*(r3,r4) (75) 

\ i=l 



Ti'(r3->r4,E) 
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(a) 





FIG. 3: (color online) Scheme of configurations giving some 
contribution to Z{t). (a) Open trajectory (OT) configurations 
counted in Z^{t). (b) Configuration with 7 and 7' surround- 
ing a periodic orbit (PO), contributing to Z'^{t). 



where i^-y' — {t^ +t^i)/2. 

The term containing the action difference is a rapidly 
oscillating function, so due to the energy average most 
of the contributions will cancel out. Only trajectories 
with similar actions will give some contribution, which 
imposes conditions on the possible configuration of the 
points r^. There are two possible configurations, as de- 
picted in Fig. 3, following the analysis in Ref. [25]: (a) 
open trajectory (OT) contributions (that we will denote 
by Z^{t)) where ri « ra and Y2 ~ r4, or additionally, 
in case of time reversal symmetry, ri « r4 and V2 ~ ra 
(this gives a factor of two, taking into account that in 
case of time reversal symmetry the eigenfunctions of H 
can be constructed to be real), (b) periodic orbit (PO) 
contributions (Z^(t)), with ri k, t2 and ra « r4 and both 
trajectories surrounding a periodic orbit. 



As before, we can calculate the 211 contribution to 
Z^{t) for K = 2. The double sum is replaced by the 
sum rule and an integral counting the encounters along 
7. The classical survival probability is modified again by 
a factor e*'""/'^''. We assume that the stability amplitudes 
of the two trajectories are the same, so the calculation of 
the integral over and Pi can be performed as for the 
diagonal approximation. Then, 



^1,211 



duds e^""w2ii(M, s, i)e-(*-*<=-)/^^ 



2e-*/^^ ( - 2— ) . (78) 



As shown before in the semiclassical evaluation of double 
sums over OT's connecting points inside a system, 'one- 
leg-loop' (111) diagrams have to be considered. The result 
for the integrals in this case is 



Z 



W=4 



t 

— ( 

tH 



(79) 



cancelling the linear contribution in Eq. (78). 

We note that this contribution can be written as 
Z^{t) ~ 2p{t), where p{t) is the mean survival probabil- 
ity of the state </'(r), i.e. p{t) = c?r|0(r,t)p. Here the 
area of integration A entering in the decay corresponds 
to the area confined by the binding potential. 

Higher order corrections can be calculated as in 
Sec. IV, and we can simply write the OT contribution 



Open trajectory contributions 

Let us consider the contributions of OT's. For this 
purpose we expand the contributions from trajectories 7 
and 7' along trajectories 7 and 7' connecting q = (ri + 
r3)/2 and Q = (ra 4- r4)/2. Thus 



Z' = np{t) , 



(80) 



where p{t) is given by Eq. (53) for the unitary case and 
by Eq. (54) for the orthogonal case. 

Periodic orbit contributions 



Z\t) 



xS{t~t 



Re ( / dQdq'dqdq' 



\ 7,7'(q-Q,B) 

^^y)A(q + q72,Q + Q72) 
xA*(q-q72,Q-Q72)e*^^-'\, (76) 



where A^-^^, ^S^- Sy - (q' • p°y - Q' • p: 



,,), and p° 



and p:^^, are the averaged initial and final momentum 
of the two trajectories, respectively. Furthermore, q' — 
ri — r2, Q' — r2 — r4 and k — I {or k — 2) in the absence 
(or presence) of time reversal symmetry. 

The diagonal approximation corresponds to 7 = 7'. To 
evaluate these terms we invoke the sum rule from Ref. 
[34], which allows us to write the integrals in Eq. (76) 



as e-*/-<^(s) \ JdqJdpAw{q,p)S{E 
in view of Eq. (71) then gives 



i7(q, p)) , which 



-t/Td 



Let us now consider the contributions of diagrams such 
as Fig. 3b. We first calculate the contribution of periodic 
orbits to the cross section, following a similar procedure 
as for deriving the semiclassical trace formula, namely by 
employing the semiclassical Green function in the defi- 
nition of a and expanding the actions around periodic 
orbits as shown in Ref. [26] 



PO 



(E) 



1 



dtAwiqj, Pi), (81) 



where the sum is over trapped periodic orbits j, Sj{E) = 
p • dq is the action integral along the periodic orbit 
and T, 
bit 



pj refers to the period of the primitive periodic or- 
f). ^ g-^'^J'^/WlTrMj 



2 1 is the stability ampli- 
tude of the PO together with the Maslov index i>j, and 
Mj is the monodromy matrix describing the lineariza- 
(77) tion around the PO. Almost all the long trajectories are 
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equally distributed in phase space if the system is ergodic. 
Therefore we approximate the time integral by a corre- 
sponding phase space average, i.e. J^" dt Aw {cli, Pi) ~ 
Tpj / drdpAw{r,p)S{E - H{r,p))/n{E), and obtain 

^PO^j^^^^M^Reyr^.D^ei^^^^K (82) 

tH . 

J 

We recognize here the form of the oscillatory part of the 
density of states. After substituting one finds that the 
contribution of periodic orbits to the cross-section form 
factor Z'^{t) corresponds to the spectral form factor of 
the open system. Substituting Eq. (82) in Eq. (70) we 
have 

(83) 

where Tjy = {Tj + Ty)/2. The expression given in Eq. 
(83) has been calculated as an expansion in t/tH in Ref. 
[30] up to 8th order for the unitary case and up to 7th 
order for the orthogonal case. In this context, Ill's do not 
play a role, since both stretches must have a minimum 
time in order to surround a PO. 

Summing up the semiclassical contributions to Z^{t), 
the decay rate (80), and Z'^{t), i.e. the spectral form 
factor Kopen{t) of the open system, we can in general 
write 



Z{t) = Kopen(i) + Kp{t) . 



(84) 



Eq. (84) is consistent with the result presented in Ref. 
[39] for t <^ tn, obtained by invoking supersymmetry 
techniques. 

For the orthogonal case this is, up to 7th order in t/tn, 



7GOE 



it): 



2 + 2^ + {N-2)L. 



(85) 



+ 



N 



+ 2 



+ 



H 



~\2 



5N 

X+3 



8\ 



''H 



19Af2 53iV 



+ 



60 

/41iV3 
V 360 

583N^ 



15 



4 



iV2 

4 



lOliV 



2520 



15 

1037V2 

60 



32 



12>2AN 32 

105 ^ y 



with N = tn/Td- For the unitary case the result reads, 
up to 8th order in t/tu. 



7GUE 



t 



+ 





N\ ^6 / 


90 ~ 








1920 


720 ^ 224 



H 

180 ^ 20 



H 



N 
28 



(86) 



Returning to the auto-correlation function by taking 
the inverse Fourier transform, we obtain for the GOE 
case 



+ 



1 1 



1 1-r^ 



Nl + T^ N'^{1 + r2)2 

{N-2) i-3r2 
iV3 (i-i-r2)3+- 



(87) 



where F = wr^. The first contribution corresponds to 
the well known Lorentzian shaped autocorrelation func- 
tion in the regime of Ericson fluctuations, first studied 
by Ericson in the context of nuclear cross-sections in 
the continuum region [40] (also experimentally observed 
[41]), and later for systems with few degrees of freedom, 
for which the corresponding classical scattering reflects 
irregular dynamics ('chaotic scattering') [42]. In the con- 
text of atomic photoionization, the Lorentzian behavior 
has been numerically [43, 44] and experimentally [3] stud- 
ied. The first and the second term in Eq. (87) have been 
derived in [25, 26], while the third term (partly of same 
order as the second one) and higher order quan- 

tum corrections to C{T) can be semiclassically assigned 
to off-diagonal loop contributions. 

For the unitary case the auto-correlation function reads 



+ 



1 1 



+ ■ 



1 1 



Ni + v^ m (i-hr2)2 

(TV -4) (1 - lor- + sr^) 



(i + r2)5 



.(88) 



Eqs. (87) and (88) are consistent with RMT results for 
indirect processes performed in Ref. [38] and with their 
expansion in powers of t/tH conjectured in Ref. [45]. In 
the following section we extend our approach beyond the 
RMT limit. 



VII. EHRENFEST TIME EFFECTS IN 
PHOTOFRAGMENTATION STATISTICS 

The Ehrenfest time te [46] separates the short-time 
quantum dynamics, where quantum wave packets follow 
the corresponding classical one, from a long-time regime 
of delocalized waves, where the dynamics is dominated 
by wave interference. Effects of this additional time scale 
have been recently considered for stationary processes in- 
volving time integration, among others, in Ref. [18, 20- 
22, 47-49]. In Refs. [16, 50] it was pointed out that te- 
signatures should be even more noticeable in the time 
domain. In Ref. [12] the T^-dependence of the leading 
quantum correction to the survival probability was cal- 
culated and provided an explanation for signiflcant devi- 
ations of numerical quantum results in the semiclassical 
regime from the RMT limit. This motivates us to extend 
our study to rg-effects in the statistics of photodissocia- 
tion cross-sections. 
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FIG. 4: (color online) Sketch of the 211 for the semiclassical 
approximation with finite Ehrenfest times. 



We follow the approach introduced in Ref. [16], for 
the spectral form factor, to calculate the Ehrenfest time 
dependence of the respective leading quantum correc- 
tions. However we distinguish, as in Ref. [22], between 
the Ehrenfest time of the closed system 



(89) 



where C is the typical system size and the de Broglie 
wavelength, and the open system Ehrenfest time. 



(90) 



related to the width w of the opening (here w corresponds 
to the number of fragmentation channels times the de 
Broglie wavelength). 

Let us consider the first (off-diagonal) quantum correc- 
tion to the correlation function C(a;) coming from open 
trajectories: 



tH 



-Re 



(91) 



As pointed out in Ref. [12] the densities should be multi- 
plied by a Heaviside function ensuring that the contribu- 
tion exists only for times larger than the encounter time. 
Only trajectories that are closer than a distance w to 
themselves will have an enhanced probability of staying. 
Correlated trajectories should come closer to themselves 
than a distance in phase space related to the opening, 
i.e., we place the PSS only in the region were the stretches 
are closer than a distance w in configuration space, see 
Fig. 4. Moreover, on the right hand side of the encounter, 
the stretches should separate at least a distance C in or- 
der to close themselves. This is because the two almost 
parallel momenta at the encounter have to grow until 
they are in exactly opposite directions, which requires 
that the stretches are no longer linearizable along each 
other and therefore should be separated by a distance 
comparable to the system size. The duration of the tra- 
jectory should then be at least 2tonc + where 



*WL = X ^ ln(£/w) 



is the time it takes for the stretches to be separated by a 
distance C when they are initially separated by a distance 



w. The weight function is slightly modified by this min- 
imal time and by ensuring that the time is long enough 
in order to have such an encounter. Thus 



w^''(u,s,t) 



(< — 2(tonc + ^Wl))' 
2r2tenc 



0{t — 2<onc — 2iwL) , 

(93) 

and the classical survival probability is modified by 
gtonc/rd^ In Appendix B the evaluation of the integral 
can be found, together with the calculation for the 111 
case. The total contribution can then be written as 



C. 



l,211+lll(' 

TE \ 



(94) 



where F = wr^, N — tn/Td and 2t|; — t'^ -\- t'^. Taking 
the Fourier transform this corresponds to a dependence 
in Z^{t) as 



rl,211-(-lll 



TdtH 



2t 



(95) 



consistent with [12] for the decay. Here we see two com- 
peting effects, on one hand if the Ehrenfest time is too 
large, loops can not be formed, 9{t — 2r|;) = 0, and there 
are no quantum contributions. On the other hand, if 
the time is long enough so that the loops can occur, 
i.e. if t > 2t|;, the probability of staying is enhanced 
by a factor e^^l'^d compared to generic orbits, revealing 
the enhanced classical survival probability due to the en- 
counter. In the energy domain, the auto-correlation func- 
tion C(w), Eq. (94), shows an exponential suppression 
of quantum effects depending on the Ehrenfest time of 
the closed system, similar to the exponential suppression 
of weak localization in transport in mesoscopic systems 
[18, 20-22, 47, 48], while additionally oscillations in to 
with a period given by r|; are expected. 

Let us consider now the Ehrenfest time dependence of 
the first quantum correction to C'^ito). A calculation of 
the Ehrenfest time dependence of the spectral form fac- 
tor of closed systems was performed in Ref. [16]. We fol- 
low here a similar approach, taking into account now the 
opening of the system, and the two different Ehrenfest 
time scales. In this situation the stretches are required 
to be separated by a distance £ also on the left and right 
hand side of the encounter. Therefore the minimal time 
for the orbits is 2tenc + 4tvvL . 

The first quantum correction to the spectral form fac- 
tor results from orbits sketched in Fig. 5 [13] denoted in 
the following by (2)^. The corrected weight function is 
then given by 



W'- > [U, S, t) = — 0(t - 2tcnc - 4iwL) 



(96) 

The contribution to the autocorrelation function, after 
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FIG. 5: (color online) Sketch of a periodic orbit with a self- 
crossing for a finite te- 



shifting the time integration by 2tcnc, can be written as 



C: 



2,(2)1 



(^) 



Re 



-{l+iuJTd)t/Td 



with 



x(t-4%L)/^'^'(w,i), 

us (t + 2tc 



(97) 



1 Iff 

TTfi Jo Jo 



xe 



^(l+2ic^ri)te„c/Td 



(98) 



The integrals can be performed as before, yielding 



Re 



-HUJT' 



E 



(99) 



(i + zr)3 Tj{i + iTY 

where T = ujTd again. Taking the Fourier transform, the 
result for the spectral form factor of the open system is 



^2,(2)1 



it) 



~t./rd„T'^/Td 



d{t 



-2 



1 + - +— r • (100) 

If Td ^ oo and the system is closed, Eq. (100) is con- 
sistent with Ref. [16]. Similarly as for Eq. (95) the step 
function ensures that only trajectories longer than 2t|; 
give a contribution, which are larger than 2r|; since the 
orbits have to close themselves. For those orbits the 
contribution is enhanced by e'^^^'^'*, again showing the 
enhanced survival probability for periodic orbits with a 
self-encounter. As in Eq. (94), Eq. (99) shows that the 
quantum corrections in the cross-section autocorrelation 
function are exponentially suppressed due to the minimal 
time that self-encounters require. In the case of periodic 
orbits, the suppression is stronger (since rj, > t^). 



CONCLUSIONS AND OUTLOOK 

We have presented a detailed semiclassical analysis of 
the quantum survival probability and of photofragmen- 
tation cross-section statistics, including higher order cor- 
rections. We have demonstrated how interference con- 
tributions associated with certain trajectory pairs pro- 
vide the key to understanding and deducing quantum 



corrections to the leading classical features in chaotic de- 
cay. We have seen in the case of the survival probability 
that the initial semiclassical treatment introduced in Ref. 
[12] for localized wave packets can be extended to non- 
localized ones by assuming a local time average, which 
allows us to treat in the arising double sums of trajecto- 
ries pairs only those that are correlated. Apart from the 
standard off-diagonal contributions, it proves necessary 
to include further, so-called one-leg-loop, diagrams in or- 
der to recover unitary, expressed via the normalization 
of the wave function when the system is closed. Trajec- 
tories with multiple encounters of several stretches lead 
to higher order corrections for systems with and without 
time reversal symmetry, for which again it is necessary 
to take into account the corresponding one-leg-loops as 
well as diagrams where both the initial and final points 
are inside encounter regions (which are not the same). 
Taking into account all the different allowed structures, 
depending on the general symmetries of the problem, we 
can reproduce RMT-type results presented in Ref. [10], 
where the survival probability was calculated using su- 
persymmetry techniques. Moreover, our approach can be 
further extended to also include systems with spin-orbit 
interaction, which corresponds to the symplectic RMT 
ensemble. We have also considered mesoscopic survival 
probability fluctuations through their variance and could 
explicitly show that they are non-universal, that is, that 
the variance depends on the spatial width of the initial 
(coherent) state. 

In the second part of the paper, we presented in de- 
tail an application of this approach to a different field, 
namely photodissociation and photoionization processes. 
We considered correlations in frequency of photofragmen- 
tation cross-sections. Its Fourier transform, the corre- 
sponding photofragmentation form factor can be semi- 
classically expressed as the sum of (twice) the survival 
probability, related to open trajectories, and the spec- 
tral form factor of the open system, related to the set 
of periodic orbits that are trapped in the open system. 
We have semiclassically computed the photofragmenta- 
tion form factor to high order in t/tu and moreover con- 
sidered Ehrenfest time effects. 

According to previous numerical results [12] there are 
clear indications for the importance of Ehrcnfest-time ef- 
fects in decay processes, leading to a shift in time of the 
quantum corrections. In the context of photofragmenta- 
tion, we have shown here that quantum corrections of the 
photodissociation form factor are also distinctly shifted 
in time (with a stronger shift for periodic orbit contri- 
butions). This time shift translates into an exponential 
suppression of quantum effects in the cross-section corre- 
lator, if the Ehrenfest time is comparable to the typical 
life time of the intermediate atomic or molecular reso- 
nant states in the fragmentation process. Our semiclas- 
sical results also predict a frequency modulation of the 
correlator with period given by the Ehrenfest time. 
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The semiclassical approach developed here to treat de- 
cay processes can be extended to address other quanti- 
ties where semiclassics so far was Umited by the diago- 
nal approximation. One example is the problem of the 
Loschmidt echo or fidelity, respectively, where a semiclas- 
sical treatment along similar lines as the one presented 
here allows one to calculate quantum corrections to the 
fidelity decay [51]. 

The present approach is still limited to times smaller 
than the Heisenberg time. An extension to longer times 
beyond tn remains as a challenging open problem of 
semiclassics for open quantum systems. 
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APPENDIX A: HIGHER ORDER 
CONTRIBUTIONS TO THE DECAY RATE IN 
THE ORTHOGONAL CASE 

Table II allows us to obtain the following semiclassical 

, t_ 

corrections to p'^ = e for the orthogonal case: 
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TABLE 11; The number of trajectory pairs and the number 
linking certain encounters for systems with time reversal sym- 
metry. 
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APPENDIX B: EHRENFEST TIME 
DEPENDENCE OF THE LEADING QUANTUM 
CORRECTION TO THE CROSS-SECTION 
CORRELATION C(w) 

Substituting the expressions (92) and (93) into Eq. (91) 
and shifting the time integral by 2*enc + 2*wl we have for 
the 211 correction to the photo cross-section correlation 

C-if '(w) = -^Re / *2 e-(i+--<i)(*+2twL)/r,j2ii(^)^^ 
tn Jo 

(107) 

with 



/-(a;) 



fJlg_ gtenc/Tdg-2(l-|-iw-r<j)tenc/T<J_ 

*enc 

(108) 

Upon change of variables, x = us/c^ and a = c/u, we 
obtain 



rA /"i 

/2''(a;) = — / dxcos{rx)x 



(l+2»wrd)/(ATa) 



(109) 



where r = (? /h, and the integral over a has already been 

performed. We compute the remaining integral by par- 
tial integration, neglecting highly oscillating terms in the 
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limit ^ ^ while keeping T^/rd (Eq. (89)) and T^/rd 
(Eq. (90)) finite. We find 



(1 + 2iuJTd) ^o,ir, -2{l+i,,r,)r%/r,. 



2Td 



(110) 



with Ehrenfest time t% = A-Hn(cVfi). Then Eq. (107) 

gives 



2t, 



C^f '(c) = - ( ^ ) e-^Re 



1 + 2iujTd _^_2i^^ 



(1 + iUJTd)^ 



(111) 

Here we used that from the definitions of r|. and ^wl 
follows + 2twL = T'j^ , and we introduced r|. = {t^ + 

A corresponding calculation can be performed for the 

111 case, where 

tn Jo 

(112) 

with 



du 



f ds f 
Jo 



X-Hn{c/\s\) 



Xgtenc/T<ig-2(l + ia;T<j)te„c/Td 



dt' (113) 



and <onc = t'+X~^ ln(c/|u|). With the change of variables 
X = usjc? and a = c/u and t" = t' + ln(c/|w|), we 
obtain 



^rTd 



7r(l + 2iwTd) Jo 



/ dxcos{rx)x^^^~^'^^'^^'^^ 
Jo 



(114) 

The integral can be evaluated as before, neglecting 
higher-order terms. Together with Cl'^^^{u)) from 
Eq. (Ill) we find 



Cifi+i"(^)=4^e-^i?e 



(1 - iujTd)^ ^-2iu,Tl 



(1 + {UTdYf 



which corresponds to Eq. (94). 
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